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We show that dressing polar molecules with a far-off-resonant optical field leads to new types of 
intermolecular potentials, which undergo a crossover from the inverse-power to oscillating behavior 
depending on the intermolecular distance, and whose parameters can be tuned by varying the 
laser intensity and wavelength. We present analytic expressions for the potential energy surfaces, 
thereby providing direct access to the parameters of an optical field required to design intermolecular 
interactions experimentally. 
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An intense far-off-resonant laser field induces a 
retarded dipole-dipole interaction between atoms or 
molecules, which is of long range character and falls off 
as 1/r, 1/r 2 , or 1/r 3 , depending on the separation r and 
the light wavevector [TJ |2]. This interaction is highly 
controllable and brings about peculiar effects in atomic 
Bose condensates, such as "gravitational self-binding," 
rotons, and density modulations leading to a supersolid- 
like behavior [3 -5 . In contrast to atoms, molecules pos- 
sess anisotropic polarizability and rotational structure, 
which renders laser-induced interactions more complex, 
and thereby offers new riches in the few- and many-body 
physics of field-dressed ultracold gases. 

In this contribution we study the interaction between 
two polar and polarizable molecules in a far-off-resonant 
laser field. We demonstrate that an optical field gives rise 
to new types of intermolecular potentials, which exhibit a 
crossover from the inverse-power decay at short distances 
to an oscillating long-range behavior, and whose param- 
eters can be altered by tuning the laser intensity and fre- 
quency. Furthermore, for a wide range of field intensities 
the problem can be described by an exactly solvable two- 
level model, which leads to simple analytic expressions 
for the effective potential energy surfaces. The theory 
is exemplified by optically-induced interactions between 
40 K 87 Rb molecules, widely employed in experiments with 
ultracold polar gases [6]. 

We consider two identical diatomic molecules, 1 and 2, 
with a dipole moment d and polarizability components, 
an and a±, parallel and perpendicular to the molecular 
axis. In a far-off-resonant radiative field of intensity /, 
the rotational levels of each molecule undergo a dynamic 
Stark shift, as given by the Hamiltonian pQ 

Hia = B^-^e 3 e* ia ^ 2 \k), (1) 

with B the rotational constant, and a l j\ h (k) the dynamic 
polarizability tensor in the laboratory frame. We as- 

* mikhail.lemeshko@gmail.com 



sume a laser beam propagating along the positive Y di- 
rection with the wavevector k = /cY, and linear polar- 
ization along the Z axis, e = Z. Given that the only 
nonzero polarizability components in the molecular frame 
are a zz = ay and a xx = Oy y = a±, and using B as a unit 
of energy, Hamiltonian (fiT) can be recast as: 

H h2 = J 2 - Arj(k) cos 2 (1>2) - v± (k), (2) 

where #(1,2) is the angle between the molecular axis and 
the polarization vector of the laser field. The dimension- 
less interaction parameter Arj(k) = r}\\(k) — r)±(k) with 
V\\,-\-(k) = a\\^(k)I / (2sqcB). We note that eq. Q was 
derived in Refs. [7J [8] using the semiclassical approach 
and the rotating wave approximation. All rotational lev- 
els exhibit a constant shift of 77^, given by the second 
term of eq. which will be omitted hereafter. 

The polarization vector of an optical field defines an 
axis of cylindrical symmetry, Z. The projection, M, of 
the angular momentum J on Z is then a good quantum 
number, while J is not. However, one can use the value 
of J of the field- free rotational state, Yj 5 m that 
adiabatically correlates with the hybrid state as a label, 
designated by J, so that | J, M; Arf) — >> Yj ? m for A77 — >• 0. 
For emphasis, we also label the values of J by tilde, e.g. 
with corresponding to J = 0. 

Induced-dipole interaction ([2| preserves parity, hy- 
bridizing states with even or odd J's, 

\J,M;At ] )=^4m( A v)Yjm, J + J even, (3) 
j 

and therefore aligns molecules in the laboratory frame. 
Aligned molecules possess no space-fixed dipole moment, 
in contrast to species oriented by an electrostatic field. 

We note that at the far-off-resonant wavelengths 
usually employed in alignment and trapping experi- 
ments (~1000 nm), the dynamic polarizability a^-(fc) 
approaches its static limit, a^-(0), for a number of 
molecules, e.g. CO, N2, and OCS. However, this is not 
the case for alkali dimers having low- lying excited 1 £ and 
x n states, such as KRb and RbCs. Virtual transitions to 
these states contribute to the ground-state dynamic po- 
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FIG. 1. Lowest rotational energy levels of a molecule in a 
far-off-resonant laser field, depending on the field-strength pa- 
rameter A rj. Different colors correspond to M — (black), 1 
(red), and 2 (blue). Energy is in units of B, two lowest tunnel- 
ing doublets are labeled as J, M. For A 77 > 15 the splitting of 
the lowest tunneling doublet can be accurately estimated as 
AE = |£|exp(3.6636-2VA^), with E = 2^A^-Ary-1 [H]. 



larizability, rendering it a few times larger than the static 
value 0[IO]. 

Figure [I] illustrates the dynamic Stark effect on rota- 
tional levels of a diatomic molecule. A far-off-resonant 
optical field of sufficiently large intensity leads to forma- 
tion of "tunneling doublets" - closely lying states of op- 
posite parity with the same M and A J = 1 [7 j. The dou- 
blet states can be mixed by extremely weak electrostatic 
fields, leading to strong molecular orientation in the lab- 
oratory frame [TTJ [12] . The energy gap between neigh- 
boring doublets increases with the field intensity, and is 
proportional to 2^/Arj in the strong- field limit [13 . As 
we demonstrate below, at large A 77 interaction between 
two ground-state molecules can be described within the 
lowest tunneling doublet. 

In the absence of fields, two polar molecules interact 
via the dipole-dipole interaction: 



V dd (r) 



(Sji - 3r jii), 



(4) 



where cr 1 ' 2 ) = d^ 1 ' 2 ) / d are unit dipole moment vectors 
of the molecules, r is the unit vector defining the inter- 
molecular axis, and 



d 2 



1/3 



(5) 



is introduced as a unit of length. In the field-free case, 
interaction between two ground-state polar molecules has 
an isotropic asymptotic behavior, V dd (r) = — l/(6r 6 ). 

Far-off-resonant laser light interacts with molecular po- 
larizability, thereby inducing oscillating dipole moments 
on each of the two molecules. Retarded interaction be- 
tween these instantaneous dipoles leads to an additional 
term in the intermolecular potential [TJ [2j [14] , 



V aa (k,r) 



•Q^f' 1 {k)Vji{k,Y)a\^ a \k)e n cos(kr), 



(6) 



with Vji the retarded dipole-dipole interaction tensor, 



Vji(k,r) = -g [(Sji — 3fjr/)(cos kr + kr sin kr) 

— (Sji — TjTi)k 2 r 2 cos kr] (7) 

For the laser light linearly polarized along the Z- 
axis, optically-induced dipole-dipole interaction ([6| can 
be rewritten in dimensionless form, 



~ l f' 2 (k) cos(kr), 



(8) 

where energy is measured in units of £?, distance in units 
of ro, and k in units of r^ 1 , and o^- = a^/Aa with 
Aa = a || — a± is the reduced polarizability tensor. The 
dimensionless parameter, 



d 2 



2Aa(fc)B' 



(9) 



characterizes the relative strength of the permanent- 
dipole and induced-dipole interactions and is on the or- 
der of 10 2 — 10 3 for polar alkali dimers. We note that in 
eq. (|6| we neglected a "static" term due to the coupling 
of the dipole moment of one molecule with the hyperpo- 
larizability of another via the optical field p~4j [15] . This 
interaction is independent of k and only becomes com- 
parable to the dipole-dipole potential Q at much larger 
intensities (/ > 10 14 W/cm 2 ). 

In the Born-Oppenheimer approximation, effective in- 
teraction potentials Kff(r) are obtained by diagonalizing 
the Hamiltonian, 



H = H 1 +H 2 + V dd + V aa , 



(10) 



for fixed intermolecular separations r = (r, 0,0). In di- 
mensionless units, the dipole-dipole potential Q is on 
the order of unity, while the strength of the optically- 
induced interaction ([8| is given by Ar?/^. Therefore, 
for the field strength parameter satisfying the inequal- 
ity 1 <C <C both interaction terms are much 
smaller than the energy gap between neighboring tunnel- 
ing doublets. Hence, in the basis of field-dressed states, 
I J1M1, J2M2), the interaction between two ground-state 
molecules can be treated within the lowest tunneling dou- 
blet, |00,00)-|10, 10). Given that the optically-induced 
potential (|8| and the dipole-dipole interaction Q mix 
only states of the same and the opposite parity, respec- 
tively p~6[ [17] , the Hamiltonian matrix takes the form: 



rel 



U dd U^ a + 2AE 



(11) 



where AE(Ar]) is the splitting between the tunneling- 
doublet levels |0,0) and |1,0), cf. Fig. [l] and the matrix 
elements are given by: 



tt / n 1 _ 3 cos 2 6> A s 
U dd (r) = = G(Ar/), 



(12) 
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FIG. 2. Short-range behavior of optically-induced KRb-KRb 
potentials in the XZ plane (0 = 0), for different values of Arj. 
Left and right columns correspond, respectively, to the ground 
|00, 00) state, and excited 1 10, 10) state potentials, as given by 
eq. |l6t . The dependence of the short-range potentials on cj> is 
negligible. Potential energy is in units of B, with V e s(r — >■ oo) 
chosen as zero. The laser beam propagates along the Y axis, 
k || Y, with the polarization e || Z. See text. 



TT j Arj(k) cos (kr) ^ N 

UL(Kr) = * J ; 3 ^ jiA^a/Aa) 

^a(kr) + [a(kr) - 36(fcr)] ^|^^) j , 

(13) 

with a = (an + 2a±)/3 the average molecular polariz- 
ability, a(x) = x 2 cosx, and b(x) = (cosx + x sinx). 

The factors G(Arj) and Kj(Ar],a/Aa) have an ana- 
lytic representation in terms of the coefficients Cjq°(A?7) 
of eq. ^ [17 . They are on the order of unity and can be 



analytically estimated in the strong- field limit, At; — > oo, 
with good accuracy: 



G(Ajy) = [l - At?- 1 / 4 ^ (iAr?- 1/4 ) 



K U0A (A V ) = y| [|i - A,-V 4 F (A," 1 / 4 ) 



(14) 



(15) 



where F(x) = exp(— x 2 ) J* ex.p(y 2 )dy is the Dawson in- 
tegral [18]. 

In such a way, the effective potentials are given by 
solutions of ( 11 ): 



V eS (k,r) = AE- 



C/L(k,r) + C/L(k,r) 



4^ d 2 d(r) + [U~UKr) - C&(k,r) + 2A25 



, (16) 



with the minus sign corresponding to the ground 1 00, 00) 
state effective potential. 

Figure [2] shows the short-range behavior of effective po- 
tentials ( |16| , induced between two 40 K 87 Rb molecules by 
a laser field of wavelength A = 1090 nm (k = 0.021/ro), 
corresponding to the following values of parameters: £ = 
96.5, a/Aa = 0.62, and r = 36.1 A EDI 021 ED] • 

The behavior of the effective potentials is dictated by 
the interplay between static dipole-dipole and optically- 
induced dipole-dipole interactions, eqs. Q and At 
small distances, kr <C 1, the dominant contribution to 



eq. (16) reads: 



V eS {kr < 1) 



II - 3 cos 2 



3 At] 



s(0) ± G(Arj) 



_ (17) 
where s(0) = sgn[l — 3 cos 2 6]. Eq. (17) reveals the pres- 
ence of a critical value, Ar] c = £y // 2/3G(A?7), that deter- 
mines the sign of the short-range potential (Ar] c ~ 65 for 
KRb). For A 77 < At7 c , eq. (17) is governed by the second 



term in the square brackets, and the potentials are always 
attractive in the ground state and repulsive in the excited 
state for any angle except for = ± arccos(±3 -1 / 2 ), 
where V e ^(kr <C 1) vanishes identically, cf. Fig [5] (a), (b). 
This behavior is qualitatively different from the dipole- 
dipole interaction between two polar molecules oriented 
along the Z axis, VjJ = (1 — 3cos 2 #)/r 3 , whose sign 
alternates in dependence on 0. On the other hand, for 
At? > A?7 C , the sign of V e s becomes angle-dependent, due 
to the interplay between the terms in the square brack- 



ets of eq. (17), resulting in the behavior similar to : 
the potential is attractive at — 0,7r and repulsive at 
= tt/2, cf. Fig. [2] (c), (d). Both below and above 
Af] c the inverse-power decay rate of V e ^ can be tuned by 
changing A 77. 

At large distances, kr ^> 1, the optically-induced po- 
tential (13) is proportional to 1/r and hence dominates 



4 




eqs. (13) and (18). Although in the case of kr <C 1 



FIG. 3. Long-range behavior of the ground |00, 00) state 
potentials in the XZ (0 = 0) and YZ (<f> = tt/2) planes, for 
Ar] = 100. Long-range behavior is similar for the excited 
1 10, 10) state, the m agn itude of the oscillations scales with 
A77, as given by eq. (18). Potential energy is in units of B, 
with V e s(r — >> 00) chosen as zero. The laser beam propagates 
along the Y axis, k || Y, with the polarization e || Z. See 
text. 



over the dipole-dipole interaction (12), resulting in an 
asymptotic behavior given by: 



V eS {kr > 1) 



3 cos(kr)^os(/cr) ^ ^ (lg) 



which manifests itself in decaying oscillations, as shown 
in Fig. [3] for the case of At? = 100. We note that the 1/r- 
like interaction (18) is of longer range than both the van 



der Waals and potentials, which may lead to intrigu- 
ing scattering properties of optically-dressed molecules. 



In general, the potentials of eq. (16) depend on the 
azimuthal angle (/>, as given by the cos(kr) term of 



the ^-dependence is negligible, the long-range behavior 
is strongly anisotropic in </>, cf. Fig. [3j As follows from 
eq. ( 18 ), the magnitude and phase of the long-range oscil- 



lations scale with At? and k respectively, and are similar 
for the ground and excited state. 

As one can see from Fig. [3j optically-induced po- 
tentials exhibit concentric minima. If these potential 
wells are deep enough, they will support long-range 
bound states, whose properties are completely deter- 
mined by the optical field, molecular dipole moments 
and polarizabilities, and are independent of the details 
of the intermolecular potential, in a similar way to 
the electrostatically-induced bound states predicted by 
Avdeenkov and Bohn [21 j. Looking into the properties 
of these states represents a challenging theoretical and 
computational problem. 

In summary, we undertook a study of intermolecular 
interactions in the presence of an intense far-off-resonant 
optical field, and provided simple analytic expressions 
for the resulting potential energy surfaces. Optically- 
induced potentials are highly controllable and are signifi- 
cantly different from the dipole-dipole interaction taking 
place between oriented polar molecules [22]. This could 
open a way to novel quantum phases of laser-dressed ul- 
tracold polar gases and new methods to control molecular 
collisions in the ultracold regime. 
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